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In the scaling regime, amplitudes for backward virtual Compton scattering and hadron annihila- 
tion into two photons depend on structure functions that describe the hadron-photon transition at 
the partonic level. We construct simple analytical models for the vector and axial- vector pion-photon 
transition distribution amplitudes based on double distributions. Via a sum rule, these models are 
l/^ . tuned to reproduce the pion-photon transition form factors. To obtain reasonable estimates of the 

' distributions at a low scale, we modify the models to saturate the positivity bounds using empir- 

\ ically parametrized parton distributions. We also determine a model-independent contribution to 

. the axial-vector transition distribution using chiral perturbation theory. 
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I. INTRODUCTION 

<^ . 

K*" . The detailed study of off-diagonal exclusive processes in perturbative QCD has lead to a wide variety of 

theoretical and experimental work. QCD factorization theorems for such processes enable the physical matrix elements 
to be written as convolutions of perturbatively calculable hard scattering amplitudes with what are now commonly 
QQ ' referred to as generalized parton distributions. Typically these distributions arise from hadronic matrix elements that 
] are off-diagonal in momentum space and contain quark or gluon operators separated by a light-like distance. Parton 
lO ' distributions, which arise in inclusive deep-inelastic scattering, involve diagonal matrix elements of light-like separated 
I operators. Elastic form factors, on the other hand, are non-diagonal matrix elements of local operators. The physics 
(-H , encompassed by generalized distributions is thus that of both inclusive and exclusive processes at the partonic level. 
O i' There are a number of insightful reviews on the subject from a variety of perspectives, see |^ Isll^ llollllj| . 
' , In the case of ordinary parton distributions, which we shall generically denote q{x), a wealth of phenomenology has 
resulted from simple analytic parameterizations. For example, one might assume at some low input scale /i the form 



q{x) = Ax-'il-xf, 

and constrain the parameters A, a, and b by particle number sum rules and fits to experimental data at higher 
scales. Typically there is considerable information buried in these few parameters, since one also generally assumes 



^ . the vanishing of sea quark distributions at the input scale, and these distributions are then radiatively generated by 
■ ■ ' perturbative evolution. Nonetheless reliable phenomenology has emerged from such Ansatze, and global analysis of 
structure function data has allowed for better tuned input parameterizations. 

For generalized parton distributions, one has functions which we shall generically write as F{x^S^,t). In addition 
to the longitudinal momentum fraction x, these functions depend on the four-momentum transfer squared, and 
the longitudinal momentum transfer, ^. Lorentz covariance imposes rather stringent constraints on the x-moments 
of F{x,£,,t), namely the n'^ moment must be a polynomial in ^ of at most n**^ degree. The coefficients of these 
polynomials are in general form factors that depend on t. For example, the zeroth moment is some form factor F{t). 
On the other hand, when the four-momentum transfer goes to zero, we have the reduction relation F{x, 0, 0) = q{x), 
which relates the forward limit of the generalized distribution to the ordinary parton distribution. One might be 
inclined to choose an input parametrization at the scale fi of the form 

F{x,C,t)^qix)Fit), 

and use existing parton phenomenology and form factor data to fix the input distribution. Without ^ dependence, the 
polynomiality property of the a;-moments it trivial. Further the reduction relation in the forward limit is guaranteed 
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by the form factor's normalization. Some problems with the resulting phenomenology are clear: one relies completely 
on perturbative evolution for the distribution to acquire ^ dependence, whereas one knows that this dependence 
arises at any scale from the interference between light-cone wave- functions of differing longitudinal momenta [T^ IT^ . 
Furthermore, there is no transparent way to improve the input parameterizations consistent with the constraints, i.e. 
once the parton distribution and form factor are known, the input is fixed. 

For generalized parton distribution phenomenology one needs a different way to construct and constrain input 
parameterizations beyond the factorized form above. In an alternate approach, one formulates the generalized distri- 
butions in terms of underlying double distributions 14, 15). The generalized parton distributions are then obtained 
as projections of these Lorentz invariant double distributions, and this feature maintains the polynomiality prop- 
erties of generalized parton distribution moments. Consequently the double distribution formulation is attractive 
from the perspective of generating input parameterizations. Calculations of model double distributions have been 
pursued flGl IT^ llsL [lp.T20ll. There is a drawback to using models based on double distributions: the positivity 
bounds [lM^[Bul2al23n^ become obscure [2^ . 

In this work we use the formalism of double distributions to provide estimates of pion-photon transition distribution 
amplitudes |2(i| |. Unlike generalized parton distributions, the off-diagonality of these transition distributions is in 
particle state, not just momentum. But dissimilar to other non-diagonal distributions considered previously p7l 
HM J the i-channel exchange involved in transition distribution amplitudes carries away nearly all the quantum 
numbers of the initial state. These distributions enter, e.g., in the pion annihilation process, 'k'^tt~ —^ 7*7, and 
backward virtual Compton scattering, j*Tr~^ —^ 7r"'"7, in the kinematic region where the real photon is nearly coUinear 
to the initial 7r+. This region is analogous to the scaling limit of virtual Compton scattering, and QCD factorization 
applies. The amplitudes for the annihilation process and backward Compton scattering factorize into a convolution 
of a hard scattering kernel with a soft matrix element, the latter describes the hadron to photon transition [26l |. 
Analogous proton to pion transition distributions amplitudes appear in the scaling limit of pN 7*7r |30f . 

This paper has the following organization. First in Sec^Jwe define the pion-photon transition distributions of [2^ in 
terms of double distributions. Next in Sec. lIIII we present a simple quark model and use it to calculate transition double 
distributions. While the model can reproduce the pion-photon transition form factors, the partonic content must be 
modified to give reasonable phenomenology for processes at large momentum transfer. In Sec. IIVI we investigate 
the positivity bounds for pion-photon transition distributions. Saturating these bounds by including empirically 
parametrized parton distributions in a factorization approximation results in phenomenologically reasonable input 
distributions satisfying known constraints. A conclusion ends the paper fSec. Iv|. 



II. DEFINITIONS 



To begin, we build up the vector and axial-vector pion-photon transition distributions by considering matrix el- 
ements of towers of twist-two operators. This approach leads us to define the transition distribution amplitudes in 
terms of double distributions. Consequently Lorentz covariance will be manifest, and this formulation will ultimately 
enable us to obtain model estimates that satisfy known constraints. 

A. Vector Operators 

The pion-to-photon matrix elements of vector twist-two operators can be decomposed in a fully Lorentz covariant 
fashion in terms of various twist-two form factors Wnk (t) , namely 

(j{P')\i^{0)T--f^fiD^'' ■ ■ ■ iB^">V(0)k+(F)) 

.--<P,A,.-l-g;g^V,.„.WP»..^P--.(--) .^.(--) . (1) 

where the action of ^ • • • ^ on Lorentz indices produces the symmetric, traceless part of the tensor, — 132 MeV is the 
pion decay constant , P is defined to be the average momentum between the initial and final states, P'^ — ^{P' + P)'^, 
and A is the momentum transfer, A'^ = (P' — P)^, with t = A^. The field ip is an iso-doublet of the up and down 
quark fields, and are the usual isospin raising and lowering operators. T-invariance restricts the form factors to 
be real, but places no further restriction on the index k in the sums. 

Ordinarily ambiguity is present in such decompositions of twist- two matrix elements [sil l32l Is^ l33 | . This situation 
stems from the generality of writing down form factors of currents that are not conserved, and leads to the freedom 
to define various distributions with the same physical content. Here, however, there is only one possible grouping of 
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Lorentz structures in the above decomposition due to parity. Nonetheless the twist-two currents are not conserved, 
and there is no ambiguity in the above decomposition into form factors Wnk{t)-^ We now define a double distribution 
W{P, a; t) that is the generating function of the twist-two form factors 

/■I /-i-l/?! 

Wnk{t)= dp da(3"~''a''W{f3,a;t). (2) 

J-l "'-1 + 1/3 

There is no a-symmetry of this double distribution because the values of k are unrestricted in Eq. 

With the operator product expansion, we can relate the moments in Eq. ^ to matrix elements of bi-local operators. 
By construction, the double distributions appear in the decomposition of the light-like separated quark bilinear 
operator 



(7(P')l^(-^/2)r-^(z/2)|7r+(P)) = ^z^^P^A.e'^^'"^ j dp j ' ' da e'^^^-^+^'^^-^/^W{p,a;t), (3) 



1 /.1-I/3I 
1 + 1/31 



where is a light-like vector. 

The light-cone correlation function is obtained by Fourier transforming with respect to the light-cone separation 

Mv{x,^,t) = ^e^-^^^-{^{P'M{^z-/2)T-^+i:{z-/2)\n+{P)). (4) 

Above, the variable ^ is the skewness parameter defined by ^ = — A+/(2P+). Unlike the case of generalized parton 
distributions, we cannot assume without loss of generality that ^ > because time-reversal does not relate the 
negative ^ distributions to those at positive The correlation function in Eq. can be written in terms of the 
vector pion-photon transition distribution amplitude V{x,^,t) 

Mv{x,U) = -=— e:P,A,£+^^-y(x,e,i). (5) 

Unlike the double distribution, the transition amplitude is a quantity that enters directly into convolutions that 
describe physical processes. Inserting the decomposition Eq. Q into the correlator in Eq. (01, we can express the 
transition amplitude as a projection of the double distribution 

/•I ri-m 

V{x,^,t)^ dp daSix- P-^a)W{P,a;t), (6) 

J-l J-l+\f3\ 

from which we interpret the ^-dependence of V{x, ^, t) as arising from different slices of the Lorentz invariant double 
distribution. 

The transition distribution satisfies an important sum rule constraint. As with the case of generalized parton 
distributions, integrating A4v{x,£,,t) over x produces the local matrix element. Thus we can relate the vector 
transition distribution to the vector form factor Fy accessible in the weak decay of the pion. Furthermore matrix 
elements of the flavor diagonal vector operators define transition distributions V" and V^, and these are related to 
the 7r° — 7 electromagnetic transition form factor, via j2^ 

J dx[QuV''{x,^,t)^QdV''{x,^,t)] = V2UF^^{t). (7) 

In our model, isospin algebra will enable us to utilize this sum rule constraint. 

B. Axial- Vector Operators 

The analysis of matrix elements of the axial- vector twist-two operators proceeds quite similarly. There are, however, 
additional contributions from so-called internal Bremsstrahlung processes that must be separated out from the matrix 



^ Electromagnetic gauge invariance is maintained for each twist- two operator as can be verified by replacing e* with P^. 
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FIG. 1: Tree-level contributions to axial-vector matrix elements in chiral perturbation theory. The pion corresponds to the 
double line, the photon to the wiggly line, while the axial current is represented by the dotted line. In the local limit, the cross 
represents the axial current and the depicted contributions are proportional to These are separated out from the structure 
dependent axial form factor. When the cross represents the axial- vector twist-two operators, the same separation can be made, 
but it involves terms proportional to /,r times moments of the pion light-cone distribution amplitude. 




elements. These contributions can be determined from chiral perturbation theory in a model independent way, while 
the structure dependent terms lead to the axial transition distribution which is modeled in subsequent sections. 

To discuss the axial transition distribution, it is useful to recall the form of the transition matrix element of the 
local axial current. For the pion transition to a real photon we have, see, e.g. [ssj ] 



(7(P')lV'(0)r-7''75V'(0)|7r+(P)) = eel 



U ( 5^'' - ) + FA{t) {P"^A'' - g'^'^P' ■ A) 



(8) 



The terms in the above decomposition proportional to the pion decay constant arise at tree-level in chiral perturbation 
theory and are depicted in Figure ^ The last term, usually called the structure dependent piece, involves the axial- 
transition form factor Fa (t) . Electromagnetic gauge invariance is only present for the structure dependent term. For 
each process, the remaining Bremsstrahlung diagrams are required for gauge invariance. The leading contributions in 
chiral perturbation theory stem from the axial- vector operator A^^ 

= ^ Tr(EV-ia^S - ET^ia^Et) (9) 

and the electromagnetically gauged axial-vector operator. Here / ~ /tt at leading order, and S is an exponential 
matrix of the pion fields 

E = exp%/2i7r • t//. (10) 

There are analogous contributions to the axial-vector twist-two matrix elements and we shall separate them out in 
order to define a structure dependent axial transition distribution. To this end, let us determine the relevant tree-level 
contribution in chiral perturbation theory. 

Following [3^ Is?! , we map the axial- vector twist-two operators onto those in chiral perturbation theory. At leading 
order^ we find 

^ 4" Y Tr(Etr'^a^^^c''^l • • • - Y.T-td^^id"' ■ ■ ■ td^-^T.''), (11) 

(n) 

where the are a set of low-energy constants which in general are not constrained by chiral symmetry. They are 



moments of the pion light-cone distribution amplitude 4>tt{x) |39|,|40j, specifically 



a 



(«) 
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1 



- dx[x--^ cj,^{x). (12) 

We demonstrate this as follows ^38| . Using the operators in Eq. pi(l , we determine 

(0|^(0)r-7<'^75iS'^i •••iB^"}V(0)k+(P)) (13) 



^ At higher orders in chiral perturbation theory, one has contributions with two extra derivatives contracted, or equivalently, by the 
equations of motion, insertions of the quark mass matrix. These as well as the one-loo p co ntributions can be systematically calculated. 
For meson distribution amplitudes, the one-loop corrections have been determined in |23|. We restrict our attention, however, only to 
those terms which survive in the chiral limit. 
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Contracting with the light-like vectors z^^z^-^ ■ ■ ■ z^^, and utilizing the operator product expansion, we arrive at 

(0|V^(-z/2)r-^75V(^/2)k+(P)) =*/.P-2 f df3e<''~^)''-^-MP)- (14) 

Jo 

And thus the familiar definition '4l] of the pion distribution amplitude emerges via a Fourier transformation, namely 

'^eK-hy-P\Q\^{-z-/2)r-M{^-l'^)\^^{P))=iUM^)- (15) 

Integration over x produces the matrix element of the axial current which defines the pion decay constant. Thus we 
require the normalization 

1 

dx4>Tr{x) = 1, (16) 



which tells us further that Og^^ — 1, see Eq. (|12|) . 

Armed with the leading-order operators (^^^i- -^" jn Eq. lll|l . we can calculate the tree-level chiral contributions 
to the axial twist-two matrix elements shown in Figure ^ The first diagram, for example, can be determined by 
electromagnetically gauging the operators in Eq. (|ll|l . and then calculating the axial pion-to-photon transition matrix 
element. One finds 

The factor of {n+X) tells us that ultimately this set of diagrams will be related to the derivative of the pion distribution 
amplitude. The leading-twist contribution from such terms, however, is a gauge artifact because z • e* = in light- 
cone gauge. Thus in the decomposition of axial twist-two transition matrix elements, we should ignore structures 
proportional to e*^. The contribution at leading twist from the second diagram is not a gauge artifact, and can be 
calculated using the pion electromagnetic current and the operators Q^^^^■■■^^" jn Eq, This calculation enables 

us to write down the decomposition of axial twist-two matrix elements that is the analogue of Eq. JH)) : 

(7(P')l^(0)T~7^''75iB^^ • • •iS^">V(0)k+(P)) 

- 2/. 4"^^^A{''(-A)^^ . . . (-A)^"> 



(18) 



where we have omitted any terms involving the structure 5'^", and used e ■ P' = for the real photon. Here the 
Bnk{t) are structure dependent form factors of the axial twist-two operators. As with the vector operators, we define 
a double distribution i?(/3, a; t) to be the generating function for these form factors 

/•I /-i-l/SI 

Bnk{t)= dp dap"-''{a + pfB{P,a;t). (19) 

J-l J-l+\l3\ 
(n) 

Using this definition and the relation of the 05 to the pion distribution amplitude, we find that matrix elements of 
the bi-locally separated quark operator have the form 

(7(P')l?(-^/2)r-^75V'(^/2)|vr+(P)) = -2eU ^^'_^^f j'^ df3 e<^'^>^M0) 



Pi fi-iPi _ 
+—P'-zA-e* dp dae-"^P-'+'''^-'/^B{p,a;t). (20) 

Jtt J-l J_i + |/3 

The light-cone correlation function is obtained by Fourier transforming with respect to the light-cone separation 



MA{x,i,t) = i / ^e"^"^"(7(P')l^(-^"/2)r-7+75V'(^-/2)k+(P)), (21) 
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and can be written in terms of the pion distribution amplitude and the axial pion-photon transition distribution 
A{x,^,t) as 

MAix,U) = eUsigniOem-\x\)^^cl,^ (^^'^ + ■ e* {1 - ^ A{x, ^,t), (22) 

where 

/•I ri~\P\ 

dp daS{x- P-^a)B{l3,a;t). (23) 

•J-l "'-1 + 1/31 

The pion pole term derived above is analogous to the one encountered in the nucleon's E{x,^,t) generalized parton 
distribution |43.l43ll43|. In our case, instead of the pion-nuclcon coupling, g-^NiO), we have the pion electromagnetic 
form factor, F(0) = (5^+ — 1. 

Integration of the amplitude ^AA{x, ^, t) over the variable x produces the plus-component of the local axial current 
operator. The normalization of the pion distribution amplitude guarantees that the pion pole term in Eq. (jSJ is 
correctly reproduced. Additionally we arrive at the sum rule constraint that relates the transition distribution of 
Eq. H22|) to the axial transition form factor in Eq. © 

J ^dxA{x,^,t)^ UFA{t). (24) 



III. QUARK MODEL FOR THE PION-PHOTON TRANSITION 

In this section, we use a quark model to calculate pion-photon transition double distributions. Our main goal is to 
obtain realistic functional dependence of these distributions on the momentum-transfer within a covariant framework. 
Moreover these models will serve as the core for constructing phenomenological estimates of the transition distribution 
amplitudes by incorporating realistic x-dependence. 

The model consists of a triplet of pions tt and a doublet of constituent quarks q. Their Lagrangian has the form 

^ = ^df,TV ■ a^TT - ^M'^tv ■ TV + q{ip-m)q~igTV -q-f^rq, (25) 

where we keep the quarks degenerate and omit color indices from the quarks, as these will produce only trivial color 
factors below. To consider electromagnetic interactions, we must add the photon's coupling to each of the particles 
in the above Lagrangian. We note that unlike the chirally invariant operators considered in the previous section, 
the Lagrangian above does not respect chiral symmetry.^ In this quark model, the quark-pion coupling g is given by 
g = mj Jt^. This model was used long ago to calculate pion-photon transition form factors ASk Alt. i48< . i49t | . Additionally 
a variety of partonic applications of this quark model have been pursued, see e.g., jsflllSlj. 

To apply the quark model to pion-photon transition distributions we must come up with a prescription to handle 
the QCD matrix elements in Eqs. (^3) and (|18|) . To make any progress in calculating the double distributions, we 
must use the parton model simplification for the QCD gauge covariant derivative: — > . Further we replace 
the elementary quark fields V' with constituent quarks q. Working in the impulse approximation, we have the quark 
model contributions depicted in FigureEl and we address the determination of the vector and axial vector distributions 
separately. 



A. Vector Distributions 



For the matrix elements of vector twist-two operators, we only have contributions from the first two diagrams. 
Explicitly the first depicted diagram has the form 

_■ /Vn * I "^'^ TrP + m)75(|^-f-Hm)7"(|^ + ^-fm)r^^^-^'"] 
^65 ^"^d^" y (27r)4[fc2-m2 + ze][(fc + A)2-m2+i£][(/fc-P)2-m2-fi£]' ^ ^ 



Our ch oice of a chirally non-symmetric Lagrangian for the constituent quarks is based on practical considerations. If one uses chiral 
quarks l45l . the pion couples derivatively and the asymptotic fall off of the pion-photon transition form factor is contrary to experiment 
and perturbative QCD. Moreover, the pion-photon transition amplitudes are complicated by divergences in the chiral quark model. For 
these reasons we have chosen Eq. 1251 as our model. 
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FIG. 2: Impulse approximation to the twist-two matrix elements. Here the twist-two operators with momentum insertion are 
denoted by a cross. The initial-state pion has momentum P, while the final-state photon has momentum P' . Only the first two 
diagrams contribute to the vector twist-two matrix elements. The last two diagrams depend on the moments of the model's 
(divergent) distribution amplitude. These plus the spectator reduced contributions from the first two diagrams must be isolated 
to deduce the remaining, structure dependent axial transition distribution. 



where Nc is the number of colors. Above we have already performed the trace over isospin and color indices. The 
symmetric, traceless tensor Y^^^'"^" in Eq. (|26|l is the momentum-space transcription of the parton model operator 
Qfifii...ii.„ g^j^^ jg given by 

pMMi-M. = -f^^'{k + A/2)^i • • • (fc + A/2)^"h (27) 

Deriving the double distribution can be done directly by casting the above contribution into the form of Eq. 1^, 
and then using the definition in Eq. |j2Jl. Taking the contribution from both the first and second diagrams, we find 
the vector transition double distribution has the form 

f 1 ^ 1 ^ 

WiP,a;t)^—lm' + -(3{a + P-l)M^-[l-a'-Pi2-(3)]-\ . (28) 

Notice in the chiral limit M^^ 0, the invariant masses of the initial and final states are the same and consequently 
Qf-symmetry of the double distribution emerges. Using isospin algebra, we can calculate the distributions for flavor 
diagonal twist-two operators. We find 

W^"(/3, a; t) = VF'^(/3, a; t) = 2W{(3, a; i), (29) 

where W{f3, a; t) is the flavor non-diagonal distribution given in Eq. I|28|) . Using these relations we can make contact 
with the sum rule in Eq. (|7|), namely 



1-/3 

dfi I daW{(3,a;t). (30) 



In the limit —^ 0, we find that the quark model reproduces the normalization of the form factor FTr-yiO) = 
(2-\/27r'^/7r)~^ required by the axial anomaly in the chiral limit |5^|U0|. Notice the variable f3 is strictly positive 
because there are no explicit sea quark contributions in this model. 

At this stage, we shall fit the constituent quark mass parameter m by requiring that the model pion-photon 
transition form factor comes close to the experimentally parametrized form. In | ,55j . a simple dipole form is used to 
describe the experimental data 

F.,it) = (31) 



where the dipole mass is A = 776 Me V. We will tune the quark mass at this stage to reproduce the form factor. 
Eventually, however, we will augment the model distribution with empirical parton distributions. This will necessitate 
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FIG. 3: Comparison of form factor fits. We compare fits of model pion-photon transition form factors to the empirically 
parametrized form in Eq. 13H . Fits for the basic model V{x,^,t) derived from Eq. 1281 and the more realistic model V^{x,£^,t) 
in Eq. H45|l are shown. 



that the model parameters be readjusted. Using the pion mass M^r = 0.14 GeV, we find that the form factor at low 
momentum transfer is well reproduced for m = 0.18 GeV. The limiting and asymptotic behaviors cannot simultaneously 
be well described with the model, however, our concern is to find viable models only at low momentum transfer and 
thus we have fit accordingly. In Figure 01 we show our fit to the experimentally parametrized transition form factor. 



B. Axial Distributions 



In considering the matrix elements of the axial- vector twist-two operators, all of the diagrams in Figure[21contribute. 
The first two diagrams give contributions quite similar to that in Eq. H26|l with the insertion of an additional 75 from 
the twist-two operator. We must be careful to remove contributions that are not structure dependent, in the sense 
of the decomposition in Eq. ||SJ). Contributions from the first two diagrams where the spectator's propagator can 
be reduced (by canceling against factors in the numerator) are not structure dependent This procedure is 

complicated by the fact that moments of the pion light-cone distribution amplitude in this model are divergent. We 
implicitly use a regularization scheme to render the distribution amplitude finite, and remind the reader of this by 
showing the third diagram in the Figure. The exact form of the scheme used is irrelevant; subsequent regularized 
contributions from the fourth diagram, and those reduced contributions from the first two diagrams are completely 
discarded because they are not structure dependent. In essence, we do not care what the model says about the pion 
light-cone distribution amplitude, for we have already determined the model independent form for this contribution 
using chiral perturbation theory, cf. Eq. H22|) . 

These technicalities aside, the structure dependent axial transition distribution amplitude can be determined from 
finite parts of the first two diagrams. We find 

B{P, a; t) ^ ^"^'^^ + -f \[3{a + P- 1) A/J - [1 - - (3{2 - ^ | , (32) 

which is very similar to the vector transition double distribution. The differences result from the differing spin 
structure factors present in the interference of the valence light-cone wave-functions of the pion and photon, and 
an overall numerical factor that depends on how the charge factors from the first two diagrams combine. For the 
parameter value m — 0.18 GeV, wc find Fa{0) / Fv{Q) — 0.98 and Fa{0) = 0.026. These are consistent with earlier 
findings for this model 48] , but are more than a factor of two too large compared with experiment [s^ . 

At this point, we do not plot the vector or axial- vector transition distribution amplitudes derived in this model. The 
main goal of this section has been to determine the <-dependence of the vector and axial form factors in a covariant 
manner. We will obtain realistic estimates of transition distributions based on these models and further considerations 
below. 
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IV. POSITIVITY AND INPUT PARAMETERIZATIONS 



The vector and axial vector double distributions determined above are presumably at some very low hadronic scale 
intrinsic to the model. It is not clear how to determine this scale from the model transition distributions alone. One 
could in turn calculate the pion parton distributions in the model, and then use the evolution equations to evolve 
empirical parameterizations down to a scale where the first few moments of the model pion parton distribution agree 
and identify this scale with the model. This procedure is not unique; many models can reproduce the empirical quark 
distributions at higher scales.'* Also the use of perturbative evolution is questionable at best at low scales. While the 
evolution kernels for generalized parton distributions are known at next-to-leading order m mil 13, 113, we will 
not use the evolution approach for our simple model. Perturbative evolution cannot generate the non-perturbative 
small-x physics which our model lacks, and we believe the small-x physics is crucial for relating to data in the scaling 
regime. The relevant leading-twist amplitudes are convolutions of transition distribution amplitudes weighted by a 
hard scattering kernel that emphasizes a region where either the initial- or final-state wave-function is evaluated near 
the end-point. In fact, the imaginary part of the amplitude is directly proportional to an overlap of light-cone wave 
functions, where either the initial or final state has a; = 0. 

In order to be useful for phenomenological estimates, a more pragmatic solution amounts to augmenting the model 
with realistic parameterizations of the pion and photon quark distributions at a low scale. We shall do this by 
implanting the realistic distribution via factorization of the /3-dependence of the double distributions |16|. This 
choice, while unrealistic for the double distributions, may indeed be less problematic for the transition distribution 
amplitudes since there is ample allowance for interplay of x, ^, and t dependence. From the perspective of model 
building, we are merely attempting to satisfy the known constraints and generate an input distribution at a moderately 
low scale. One can then evolve upward to the scale relevant to make contact with data and eventually learn how to 
improve the input parametrization. 

To implant realistic partonic behavior and estimate the transition distribution amplitudes, we shall use a factoriza- 
tion Ansatz motivated by the positvity constraints. Just like generalized parton distributions, the transition distribu- 
tion amplitudes satisfy constraints arising from the norm on Hilbert space: the positivity bounds (12 L flS l l2ll i2i. |23L l23| . 
Using the method described in , we can derive the basic bounds satisfied by both the vector and axial vector am- 
plitudes. For the positively charged pion, we find 



MvA^.tt) <\ld'< [^ju-^ l^^]. (33) 



Here we have used f^{x) to denote the parton distribution of the flavor / in the state H. The vector and axial- vector 
distributions derived in the quark model above self-consistently satisfy the positivity bounds. That is, the constraints 
in Eq. H33|l are met using the quark distributions of the pion and photon determined within the model. ^ For useful 
phenomenology, however, these model parton distributions must be modified, and hopefully in such a way that the 
bound in Eq. (|33|l remains satisfied. As an inspired guess, we choose to saturate the bounds in the forward limit by 
making a factorization Ansatz of the form 

W{(3,a;t) 



W''iP,a;t)=Nv JdHP)u^^{P). (34) 



for the vector distribution and 



B{P,a-t) 



B^{(3,a;t)^NA ' ] ' -d d^^ {P)u-y(3), (35) 

SJlpdaB{P,a-{)) ^ 



* There is no guarantee, moreover, that the scale of the imphcit quark-distribution of the photon is not mismatched with the determined 
scale of the pion's quark distribution. One should thus calculate the model's quark distributions for both the pion and photon and 
determine their mutual scale by a simultaneous fit. 

^ This can be easily understood due to the simplicity of the quark model used. When one confronts the light-cone correlation functions 
Mv A directly instead of via the double distributions, one finds that the contributions for x > ^ stem from the light-cone time-ordering 
where the spectator quark is on shell. The remaining two propagators then reduce to the model light-cone wave-functions for the initial 
and final states. A simple application of the Cauchy-Schwarz inequality gives the bound in Eq. 1331 . because the quark distributions 
are transverse momentum integrals of the squares of light-cone wave-functions. 
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for the axial- vector distribution. Here Ny and Na are constants to be determined below. When {t, ^} 0, integration 
over P and a of the Ansatze yields 



A^(x,0,0) = NA\/di (x) (x), (36) 

where V^^(a;,^,i) and A^(x,0, 0) are the realistic vector and axial-vector transition distribution amplitudes formed 
from the factorized double distributions. From Eq. H36|l . one sees our attempt to input realistic light-cone wave- 
functions for the pion and photon consistent with positivity. The resulting form is reasonable provided there is little 
interference arising from the transverse momentum components of the wave- functions. We are not advocating that 
this is the case. We are only estimating the distributions to the best of our ability given reasonable x-dependence of 
parton distributions. Improving these estimates sensibly is difficult and the subject of an ongoing investigation. In 
fact for our model, the a integration in Eq. H35I) introduces a singularity proportional to 1/P in the double distribution. 
This suggests that interference between the various Fock components in the axial distribution cannot be neglected. 
In our model for the axial distribution, saturation of the positivity bounds in the forward limit is not possible. To 
tame the singularity, we shall add a factor of /? to the numerator of Eq. as the simplest way to introduce the 
requisite interference: 

B^{P,a;t) = Na ,1,^^"''^ JdHP)u-Uf3). (37) 

In the literature, there exists a number of analytic parameterizations of quark distributions that have been con- 
strained to fit data. Our preference is to choose those models which are at a reasonably low scale, so that evolving 
up to the relevant scale for comparison with experiment allows some of the crudeness in the factorization Ansatze to 
be evolved away. For the pion, in Ref. simple analytic parameterizations of the valence and sea quark distribu- 
tions are presented. As the sea quark distributions are unconstrained experimentally, they have been related, using a 
constituent quark model, to the radiatively generated parton distributions in the proton, see ^63, 64]. For our model 
estimates, we need the up distribution in the positively charged pion, 

u"^(x) = iv"^(x) +r^(2;). (38) 

This distribution is given in terms of the valence and sea distributions defined in that are given by 

z;"^(x,//2 ^0.40Gev2) ^ 1.391 x-^^-^^^l - xf -^^"^ , (39) 
q^^{x,fi^ = OAOGeV'^) = 0.4172:-" ™(1 - 2.466 4- 3.855 a;)(l - x)'''^^^. (40) 

The input scale /i^ = 0.40 GeV^ is obtained from using next-to-leading order QCD evolution. 

For the real photon, we use the analytic parameterization presented in 65]. The photon's quark distribution consists 
of a perturbative point-like contribution and a non-perturbative hadronic contribution. The former is chosen to vanish 
at the input scale = 0.40 GeV^. The hadronic piece is determined using a vector meson dominance model |6^, and 
by relating the quark distributions in vector mesons to those in the pion. For the down quark distribution of the 
photon 

d''{x) = l-Gld-\x), (41) 

where « 0.25, and d'" ~ ^(d'"^ + d'^ ). In terms of the parametrized forms for the pion parton distribution used 
above, d''" = +4^^^^]. 

The terms motivated by the positivity bound in the region ^ < a: < 1 all concern the quark distributions of the 
pion and photon and hence contribute to the transition distribution amplitudes in the kinematic region — ^ < x < 1. 
Assuming the distributions vanish outside this range is inconsistent with the partonic input. To arrive at the anti- 
quark contributions in the kinematic region — 1 < a; < ^, we make a similar set of Ansatze by considering the positivity 
bound in the region — 1 < a; < In this region, we have 



< Jrf^ ( ^r|) (^V (42) 
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FIG. 4: Vector pion-photon transition distribution amplitude. The function V'^{x,(,,t) is plotted at fixed t = — O.lOGeV^ for a 
few values of ^. Notice the distribution is only qualitatively x anti-symmetric. 



In terms of the parameterizations used above, u"^ {x) — If^ [x) and <P(x) = d?{x). For lack of more information, we 
are forced to assume the anti-quark contributions have the same t-dependence as the quark contributions, and thus 
take the anti-quark distributions to be of the form 



W{(i,a-t) 
!]_~lpdaW{l3,a-0) 



W"-{(3,a-t)=Nv -Vrf''(/3)«"^(/?), (43) 



for the vector distribution and 



PB{P,a-t) 
fl1pdaB{l3,a\Q) 



B'\P,a-t)=NA .S^'^'^y (44) 



for the axial-vector distribution. A factor of /3 has been inserted because the axial distribution becomes singular if 
one saturates the positivity bound in the forward limit. We have defined these functions for positive /3 but have used 
a barred notation to remind the reader that they are anti-quark contributions. The realistic transition distribution 
amplitudes are formed from the barred and unbarred double distributions via 

/•I /-i-zs 

V^\x,i,t) = df3 [S{x-P-^a)W^{P,a;t)-6{x + f] + ^a)W^{(3,a;t)], (45) 

Jo J-l+0 

l-l 1-1-13 _ 

A"{x,£„t) = d(3 [S{x- (3-£,a)B'^iP,a;t)-6{x + P + ^a)B"i(3,a;t)], (46) 

Jo J -1+13 

where the minus sign arises from the opposite ordering of the anti-quark creation and annihilation operators, see |12| . 

Using these parameterized parton distribution functions in the Ansatze Eqs. H34(l . (|37|l . 143|l . and (|44|l . we can 
determine the constants Ny and Na by using the experimental numbers for -FV(O) derived from CVC, and -Fa(O). 
Using Fv{0) = 0.0259, and ^^(0) = 0.0115 iHg, we find Ny = 0.33 and Na = 0.48. There is no reason to suppose 
that the constituent quark mass parameter m remains the same when augmenting the double distribution with 
parametrized quark distributions. Just as we did for the basic model, we can determine m for the realistic model by 
requiring a reasonable fit to the small momentum transfer behavior of the pion-photon transition form factor _F!n.^(i). 
In Figure 131 we have shown the form factor determined from V^{x,S^,t) for the value m = 0.20 GeV, as well as the 
earlier fit, and the empirical dipole form. We are able to match the low momentum behavior of the transition form 
factor. Having determined the model parameters for the vector and axial vector distributions in Eqs. (|45|1 . and l|46f) . 
we can now plot the estimates of transition distributions amplitudes. The vector distribution V^{x,^^t) is plotted as 
a function of x at fixed t = —0.10 GeV^ for a few values of ^ in FigureQl The same is done in FigureOfor the function 
A^{x,£,,t). One should note that in the symmetric kinematics, we have the relation 

* = — x-e ^^^^ 

Thus for fixed i, we must have 
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FIG. 5: Axial-vector pion-photon transition distribution amplitude. The function A^{x,^,t) is plotted at fixed t = — O.lOGeV^ 
for a few values of ^. The distribution is only qualitatively x anti-symmetric 



and unlike the case of deeply virtual Compton scattering, the maximal skewness is not kinematically restricted. In the 
Figures, the minimal skewness at t = — O.lOGeV^ is ^min = —0.72. We have not plotted the distributions for negative 
values of ^. In our model the differences between distributions at positive and negative values of ^ are minimal. The 
approximate ^-symmetry that emerges is due to the factorization Ansatze which evaluate the pion and photon parton 
distributions symmetrically. 



Above we have written the vector and axial-vector pion-photon transition distribution amplitudes in terms of 
transition double distributions. These double distributions were defined from the matrix elements of the relevant 
twist-two operators. For the axial-vector transition distribution, we determined the form of a model-independent 
contribution depending on the pion light-cone distribution amplitude from chiral perturbation theory. 

Using a quark model, we calculated the transition double distributions, and demonstrated that the fall-off of the 
pion-photon electromagnetic transition form factor could be well reproduced by tuning the constituent quark mass 
parameter. Motivated by the positivity bounds, we then took a factorized Ansatz for the double distributions in 
order to establish the renormalization scale of the model, and to incorporate parameterized parton distributions with 
realistic x-dependence. Model parameters were then re-tuned to produce the pion-photon electromagnetic transition 
form factor and the pion axial form factor. We then plotted the resultant input distributions which should lead to 
reasonable phenomenological estimates for pion-photon transitions in the scaling regime. It would be interesting to 
compare our estimates with other calculations. There are a number of existing realistic studies of pion structure that 
could be extended to transition distributions amplitudes, e.g. ^67, 68. .69„.70, 71j . 

Improving the factorization Ansatz employed above may ultimately be necessary to best describe the x, ^, and t 
dependence of the data. There is ample room for improvement, e.g., we assumed the same ^-dependence and absolute 
normalization for quark and anti-quark contributions. More realistic models, data from lattice QCD, and experimental 
data can help to motivate better input parameterizations. Furthermore the only difference between vector and axial 
vector distributions in our model arises from the spin structure of the valence light-cone wave- functions, whereas 
realistically the difference stems from the interference from all Fock components. One needs to address whether 
our approximation is sensible for the input distributions at a low scale. We already found suggestions of sizable 
interference between Fock components in the axial distribution. Finally approximate ^-symmetry of the distributions 
emerges due to the simplicity of our Ansatze. Better input parameterizations must address how to treat the initial 
and final states differently in a covariant fashion. 

The study of off-diagonal exclusive reactions leads to phenomenology involving generalized parton distributions. 
Unlike ordinary parton distributions, input parameterizations cannot simply be written down at will, as constraints 
on the distributions are non-trivial. Using double distributions, we have formulated a basic input parameterization 
for pion-photon transition distribution amplitudes. Contact with future data, lattice QCD data, or other models 
will help one to tune this type of input parameterization. This procedure in turn gives one a better sense of how to 
construct reasonable input parameterizations of generalized parton distributions in general. 



V. 



CONCLUSION 
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